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Abstract 

We derive the semiclassical approximation to Feynman's path inte- 
gral representation of the energy Green function of a massless particle in 
the shadow region of an ideal obstacle in a medium. The wavelength of 
the particle is assumed to be comparable to or smaller than any relevant 
length of the problem. Classical paths with extremal length partially creep 
along the obstacle and their fluctuations are subject to non-holonomic con- 
straints. If the medium is a vacuum, the asymptotic contribution from a 
single classical path of overall length L to the energy Green function at 
energy E is that of a non-relativistic particle of mass E/(? moving in the 
two-dimensional space orthogonal to the classical path for a time r = L/c. 
Dirichlet boundary conditions at the surface of the obstacle constrain the 
motion of the particle to the exterior half-space and result in an effective 
time-dependent but spatially constant force that is inversely proportional 
to the radius of curvature of the classical path. We relate the diffractive, 
classically forbidden motion in the "creeping" case to the classically al- 
lowed motion in the "whispering gallery" case by analytic continuation in 
the curvature of the classical path. The non-holonomic constraint implies 
that the surface of the obstacle becomes a zero-dimensional caustic of the 
particle's motion. We solve this problem for extremal rays with piecewise 
constant curvature and provide uniform asymptotic expressions that are 
approximately valid in the penumbra as well as in the deep shadow of a 
sphere. 

PACS: 03.65.Sq, 42.25.Fx, ILlS.Kc 

1 Introduction 

It was some three hundred years after Descarte's work on diffraction theory that 
Keller Q provided a sound mathematical foundation for the theory by examining 
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the asymptotic expansion of the energy Green function. Armed with that sohd 
footing, and using considerable physical insight, Keller and co-workers |2 El and 
others^] found a host of applications. Though the theory has had any number 
of successes, its application can sometimes be cumbersome and its validity is 
limited to the deep-shadow region. The limitation is not really a restriction 
in the asymptotic evaluation for high energies, since any point in this case is 
either in the lit or the deep shadow region of an obstacle. It does however point 
to the limitations of the approach itself that have only been overcomefQ by 
considering a different expansion in the penumbra. 

All these approaches are purely classical, based on the eikonal- (or ray-) ap- 
proximation to the wave equation. The asymptotic forms of the wave function 
and of the energy Green function thereby obtained give the diffraction pattern 
for wave number k oo, providing the leading correction to geometrical op- 
tics due to diffraction. Keller interpreted his result P geometrically in terms 
of paths, including those of certain "creeping" rays, a very physical and useful 
way of viewing the problem. We here study diffraction from the point of view 
of a semiclassical expansion. Our starting point is the Feynman path integral 
representation for the Green function of a massless scalar field in the presence 
of ideal obstacles. To simplify the presentation we assume Dirichlet bound- 
ary conditions for the field at the surface of the obstacles and often assume a 
uniform surrounding medium. We analyze and approximate the path integral 
representation of the Green function semiclassically. Note that Feynman's ap- 
proach was published in 1949|7j, but its advantages for obtaining asymptotic 
solutions - not only technically, but in the clearer physical picture it provides 
- may not have been apparent at the time the modern version of diffraction 
theory was developed. Furthermore, it seemed natural to attack a classical the- 
ory classically. In fact, Gutzwiller[S] and others have made significant advances 
in the semiclassical evaluation of the path integral. To our knowledge these 
advantages of Feynman's formulation and Gutzwiller's semiclassical evaluation 
have never been fully exploited in the case of diffraction theory. (Some related 
approaches by others will be briefly commented on below.) We here hope to 
help close this gap. 

That an asymptotic analysis of the classical wave equation for large k should 
give the same result as the semiclassical expansion is well known but may nev- 
ertheless warrant a remark. For simplicity, we assume for the moment that the 
medium is uniform. Since the connection between the local wave number k and 
the energy £^ of a massless particle is 

k = E/{hv{E)) , (1) 

where v{E) is the phase velocity in the medium, letting k be arbitrarily large is 
equivalent to letting Ti be arbitrarily small. 

Our final result is equivalent to those obtained previously in certain limits. 
However, the present development may have a number of conceptual advan- 
tages. The illuminated region, the penumbra and the umbra are all treated in a 
unified fashion. The approach furthermore is probably more accessible to most 



2 



physicists than the original one, and might provide new insights. By providing 
a different point of view of the phenomenon, it may also suggest other (inter- 
polating) approximations. Perhaps the most interesting aspect of the present 
formulation is that it provides a systematic geometrical framework that in prin- 
ciple allows one to obtain asymptotic expressions for diffraction for many sys- 
tems, although considerable numerical work may be necessary if the geometrical 
setting is sufficiently complicated; in particular, it shows that diffraction can be 
understood as a special case of quantum mechanical tunnelling in the sense that 
the corresponding motion is classically forbidden and that the semiclassical ap- 
proximation is obtained by analytic continuation in some of the parameters of 
the problem. 

We will study the semiclassical approximation, G(Y,X;£^), to the energy 
Green function 

g{X,Y;E) = (Y\g{E + ifj)\x)^(Y\{H-E-ir^)-'\x) , (2) 

the coordinate representation of the Laplace transform of the time evolution 
operator. For simplicity we will study only the propagation of a free massless 
scalar particle in a medium which, for the most part, will be taken to be ho- 
mogeneous. We require that Q vanish for |X — Y| — > oo and also whenever X 
or Y are on a smooth compact (but not necessarily connected) two-dimensional 
surface J-^. These Dirichlet boundary conditions are also the appropriate ones 
for a particular polarization of the electromagnetic field in the presence of ideal 
conductors. Note that semiclassically only photons whose polarization is normal 
to the surface can diffract in the case of an ideal conductor. 

Green functions are among the most basic concepts in physics. Thus, for 
example, the trace oi Q{E + irj) is the response function g{E), whose imaginary 
part gives the spectral density g{E), where g{E)dE is the number of states with 
energy between E and E -f dE. Unfortunately, Green functions can rarely be 
obtained exactly. An exceedingly useful result, largely due to Gutzwiller|S], is 
the observation that the semi-classical approximation to Q{Y ,X; E) and g{E) 
are usually completely determined by the solution of a corresponding classical 
problem. Unfortunately, the solution must be modified if it is to be used in the 
case of diffraction. The corresponding classical trajectories are non-holonomic 
constrained extrema of the action that are not stationary; the dependence of 
the action on the fluctuations therefore contains linear terms. A naive ap- 
plication of the usual semiclassical approximation fails. We will see that the 
semiclassical Green function G(Y, X; E) nevertheless is given by the solution 
to a corresponding classical problem in this case too. The corresponding clas- 
sical problem one has to solve in the case of diffraction is that of the motion of 
a non-relativistic particle in two dimensions under the influence of a generally 
time-dependent gravitational-like acceleration inversely proportional to R{s{t)) 
in the presence of a ceiling; R(s{t)) is the radius of curvature of the constrained 
classical trajectory connecting X and Y as a function of the arc length s{t). 

Our interest in semiclassical diffraction theory was rekindled in the course 
of a study of Casimir effects. We found that a semiclassical evaluation gives 
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the leading behavior whenever the Casimir energy diverges |H] as one length 
scale of the problem is taken to be much larger than any other. In many cases 
this asymptotic evaluation is sufficient to determine the Casimir energy exactly. 
A phenomenologically interesting example is Derijaguin's problem|lU| of the 
Casimir force between two conducting spheres of radii Ri and R2 in the limit 
where their separation d is arbitrarily small compared to i?i and i?2- This 
force had been determined using very plausible arguments |10l[TT] . but the semi- 
classical calculation provided what we believe to be the first rigorous derivation 
of this result. The Casimir force in this limiting case is proportional to R/d^, 
where R = RiR2/{Ri + R2), and diverges as R/d 00; the semiclassical 
calculation gives the correct coefficient. On the other hand, for separations large 
compared with the radius of either sphere, the Casimir-Polder forceIT21 falls off 
as 1/d^ and a priori one cannot assert that the semiclassical approximation 
must be exact in this case. To see whether this approximation qualitatively 
reproduces the 1/d^ behavior of the force in this limit, one has to take diffraction 
into account at least semiclassically. The effort to include these effects and thus 
perhaps extend the validity of the semiclassical calculation beyond the regime 
d<icR led us to the present study. The extension to d ^ ^ is under consideration. 

In Sec. II we present a heuristic argument based on Fermat's principle to de- 
rive the effective action for the semiclassical approximation to the path-integral 
representation of the energy Green function of a massless particle in the presence 
of idealized obstacles. There is no stationary classical path that can describe 
propagation into the shadow region of the obstacle. The situation turns out to 
be analogous to that of calculating the amplitude that a non-relativistic par- 
ticle move in finite time from the ceiling to the ceiling under the influence of 
a gravitational force. We show in Sec. Ill that an analytic continuation in the 
parameters solves the problem by changing it to that of a non-relativistic par- 
ticle moving in finite time from the floor to the floor under the influence of 
a gravitational-like force - the same problem encountered in the semiclassical 
description of "whispering galleries". The remainder of Sec. Ill is devoted to 
an analytical study of diffraction by a sphere. In Sec. IV we relate the descrip- 
tion by classical trajectories in a gravitational-like field to Keller's asymptotic 
evaluation of the energy Green function for this example. In particular, we also 
derive the power corrections of the asymptotic expansion within this unified 
approach. 

There is an immense literature on the subject. See, for example GutzwillerjH] 
and ReichlJ2|j and references therein. We also include a small sampling of 
recent developments on caustics [Tl]. 

2 The Action: A Heuristic Argument 

Our ultimate goal is to obtain the semiclassical approximation G(X, Y; E) to 
^(X, Y;iJ), the probability amplitude that a photon of energy E at the initial 
point X will at some time appear at Y. Since the semiclassical energy Green 
function is essentially determined by the quadratic action S"""^ for the fluctu- 
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ations - as we will discuss later - we first determine the form of this action. 
(See Eq. (HTj) below.) 



2.1 The Action for the Fluctuations to Quadratic Order 

We argue that the contribution to the energy Green function of a continuous 
path 7 from X to Y is determined by the action 

5(7, E) = Et{^, E)=E r dt\±{t)\/v{X{t),E) . (3) 
Jo 

A path 7 from X to Y in this context is defined by the three continuous functions 
on the interval Q < t < T that make up the vector X(t) and satisfy X(0) = X 
and X(r) ~ Y. The variable t G [0,T] need not be the actual time, since the 
integral in Eq. ||2J) does not depend on how the path is parameterized. Using 
the arc length s along the path 7, one may rewrite S{"f, E) as 

where is the path's total length. For the semiclassical quantization it will be 
important that Eq. in fact does not depend on the "speed" |X|, subject to 
the restriction that |X| never vanish along the path (or equivalently that the arc 
length is a monotonically increasing function of the parameter t). 2t:i>t{'-/, E) is 
the phase lag of a monochromatic wave of frequency v — E / {2t:K) along the path 
7 with the phase velocity v{s, E) related to the local index of refraction n{s, u) = 
c/v{s, E) of the medium along the path. Fermat's principle is the statement that 
the "classical" trajectory 7c of a monochromatic ray is an extremum of r. It 
is the phase velocity that determines the trajectory of a monochromatic ray 
(as is for instance evident from Snell's law). Other than when the medium is a 
vacuum, the time t(7c, E) is not the time between emission at X and absorption 
at Y of a photon of frequency v. The travel time for a photon or, classically, a 
wave-packet, is determined by the local group velocity Wg(s, E) along the path, 

, (5) 



dE \v{s,E) 

that is, the time T for a photon of energy E to travel from X to Y is 

Jo Vgis,E) dEJo %is,E) dE - ^' 

The relation Eq. ^ between the total time T and the quantity S defined 
by Eq. (PI, in conjunction with Fermat's principle, identifies S^jjE) as the 
classical action that describes the motion of a massless particle (uniquely up to 
an irrelevant constant that does not depend on the energy nor on the path 7). 
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Hamilton's principal function i?(X, Y; T) with the independent variable T is 
the Legendre transform of S, 



i?(X, Y; T) = 5(X, Y, E)-ET=E^^ 



ds 



,E) 



(7) 



E^E{T) 

and vanishes in a vacuum, as it must, because the Minkowski distance between 
the events vanishes and Hamilton's principal function is proportional to it. 

We temporarily restrict our considerations to the situation where the phase 
velocity v{l^{t),E) — v{E) is constant when X(t) is exterior to the volume 
V with surface J- and vanishes when X(i) is in the interior of V. This is an 
idealized limit of the physical situation where obstacles with a very high index 
of refraction are embedded in a homogeneous medium. The index of refraction 
in the physical case depends smoothly on the coordinates and the classical paths 
are stationary points of the action. The action defined by Eq. Q of a path that 
avoids V in going from X to Y is proportional to its length. The vanishing phase 
velocity in the interior of V leads to a non-holonomic constraint on the classical 
motion. We will see that the situation is akin to motion in a gravitational-like 
potential in the presence of a ceiling. The classical motion is still an extremum 
of the action, simply because the constraint can be viewed as a particular limit 
of the physical case. In the limit, the extremum generally is not a stationary 
point of the action. The change of the action under some small deviations from 
the extreme path can no longer be made arbitrarily small in this limit, that is, 
the functional derivative of the action at the extremum does not exist. 



X 















'V 





Fig. 1: The semiclassical expansion. The bold curve is a schematic plot of a classical 
path 7c of extremal length from X to Y that is excluded from the (smooth) volume 
V. The length of 7c is extremal but not stationary, since it "creeps" along the surface 
between the points Xo and Yq. A non-extremal path 7 in the vicinity of this classical 
path is shown as a thin line and the local coordinate frame at the point Q(t) of 7c is 
sketched. The fluctuation 7 is fully described by the coordinates Xr(t) and x±(t) of 
local coordinate frames along the whole classical path. Note that lies in the surface 
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As can be seen from Fig.l, the classical trajectory of extremal length from X 
to Y can be rather complicated even in the presence of smooth obstacles. It may 
consist not only of straight sections and reflections off points on the surface, but 
also of segments that "creep" along the surface. The classical trajectory is by 
definition continuous and of extremal length in coordinate space. It is natural to 
parameterize it by the "intrinsic" time t for which s{t) — v{E)t is the arc length 
along the classical path. In this case the instantaneous velocity Xc(t) and the 
instantaneous acceleration X it) are orthogonal vectors, since |Xc(t)| = v{E) 
does not depend on time. 

This singles out a natural local orthonormal coordinate frame at each point 
along curved sections of the classical path. The basis vectors of this local coor- 
dinate frame are 

= = ' e^W=e'IWxe'-(i) , (8) 

where the local radius of curvature of the classical path is related to the in- 
stantaneous acceleration by R{t) = w^(i?)/|X(t)|. ell(i) clearly is parallel to 
the velocity X(t) and thus tangent to the classical path, e''(t) lies in the local 
plane of motion (and is normal to the surface of the obstacle) while e^{t) is 
perpendicular to the local plane (and lies in the surface of the obstacle). On 
straight sections of the classical path a similar local orthonormal coordinate 
frame exists, but is not uniquely defined with respect to rotations about the 
classical path. 

We next introduce local coordinates that give the deviation of the path 7 
from a classical path 7c that is an extremum of the action in Eq. Let Xc(t) 
and X(t) describe a point on the classical trajectory 7c and a point on the path 
7 in a (time-independent) global coordinate frame at the intrinsic time t of the 
classical path. The coordinates x\\ (t), Xr{t) and x^{t) for the fiuctuations in the 
local basis at time t are then defined by 

X(i) =Xc(t)+a;||(t)ell(t)-f a;^(i)e'^(t)+a;i(i)e^(t) . (9) 

where X(t) is a point on 7. (See Fig. 1). It is important to note that certain 
fluctuations just reparameterize the classical path. This follows from our earlier 
observation that the speed along the path plays no role. More formally, if instead 
of the intrinsic time t one had chosen t', where 

t'^t + r,it) , (10) 

to parameterize the classical path, then to first order in 77 (i) 

Xc(t') = Xc(<) + r?(<)Xc(i) + . . . = Xc(t) + 7^it)v{E)J it) + ... (11) 

An infinitesimal longitudinal fluctuation, x\\(t) = r](t)v{E), thus describes the 
same classical path and is equivalent to the use of a slightly different "clock" 
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to parameterize the classical path. To avoid this ambiguity and have a one- 
to-one correspondence between the fluctuations and paths, one considers only 
fluctuations in the planes 

X||(t)=0. (12) 

We can do so, because the local coordinates Xr{t) and x±{t) in each plane 
along the classical path 7c completely specify 7. The geometrical meaning and 
uniqueness of this construction can be seen in Fig. 1. 

We will also need to know how the basis vectors defined by Eq. ^ rotate 
with time. Since, by Eq. ||SJ), 

ell(t) = X{t)/v{E) = -v{E)e'-{t)/R{t) (13) 

is in the direction of one of the basis vectors, the rotation of the coordinate 
frame with time is specified by only two (generally time-dependent) angular 
velocities, instead of the usual three. The basis at time t is related to the basis 
at time t + dt hy an infinitesimal orthogonal transformation Q. ell(t 4- dt) thus 
is given in terms of the elements of and of the e(t)'s of Eq. 0. Comparing 
the resulting expression with Eq. Ijl^f) determines two of the three independent 
elements of fl. Denoting the as yet to be determined independent element of 
by a{t), one obtains 

^"^^^^"Wy^''^^)' e^it)^a{t)e^t). (14) 

Explicitly differentiating e^{t) in Eq. (jS)) and taking the inner product with 
e-^{t) determines a{t) as 

«W = ^(Xc(t)-e^(t)). (15) 

a{t) vanishes if the classical trajectory lies in a plane and thus is a measure 
of its skewness. It is the rate at which the coordinate frame rotates about the 
tangent vector to the classical path. For a helix on a cylinder of radius i?cyi, 
a velocity Vz parallel to the cylinder axis and an angular velocity u, one finds 
a{t) — —LoVz/y/v'^ + {ujRcyi)^, independent of time 

We are now in a position to expand the action of Eq. ^ to quadratic order in 
the fluctuations a;|| = 0, Xr(t), xj_{t) that uniquely describe a path in the vicinity 
of the classical one. Taking the time derivative of Eq. jnj and using Eq. (|14|l to 
describe the time dependence of the basis one finds 

X(i) = [{±r{t)+a{t)x±{t))e''{t) + {x±{t)-a{t)xr{t))e^{t)] 

+v{E) + ^yHt) . (16) 

As noted earlier v{E)e''^{t) is the velocity Xc(t) along the classical trajectory. 
[ Eq. gives the same relationship, but with Xc(t) rather than X(i). There 
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is no inconsistency, since we have chosen xi{t) = 0.] Expanding the action 
of Eq. (PI to quadratic order in the fluctuations using Eq. (|16|l . one finally 
obtains 



t(7c,-B) 

dt 



E + E 



Xr{t) 



+ 



E 



2v^{E) 



{xr{t) + a{t)x^{t) f + {x^{t) - a{t)xr{t)y 



(17) 



for the action of the path 7 in semiclassical approximation, when the classical 
path is 7c. Note that the term linear in the fluctuations of the semiclassical 
action H17|) is due to the fact that the classical path, constrained to lie outside 
of V, is extremal, but not stationary. Since the classical trajectory described 
by Xr{t) = x±(t) — creeps along the surface whenever < R{t) < 00 and 
since e''(i) is normal to at that point, the coordinate Xr{t) takes only positive 
values on a creeping segment. The kinetic terms in Eq. H17() are those of a 
non-relativistic particle of mass 



ruE = E/v'^{E) 



(18) 



moving in a plane that is rotating with angular velocity a{t). To explicitly see 
this, we perform a time-dependent orthogonal transformation of the coordinates 



Xr{t) = x^{t) cosd{t) — Xj^{t) sm9{t) 
x±{t) = x^{t) sine{t) + Xj_{t) cos 9{t) 



with 



9(1) 



dt'a{t') 



The expressions in Eq. (|17|l are then 

Xr{t) + a{t)x±{t) 
x±{t) — a{t)xr{t) 



x^{t) cos 6{t) — Xj_{t) sm9{t) 
x^{t) sin 9{t) + xj_{t) cos 9{t) 



(19) 
(20) 

(21) 



Written in terms of x^ (t) and x^ (t) the semiclassical action of Eq. ()17|l becomes 



(22) 



where the two-dimensional Lagrangian, 



E,t) 



E+^{xl{t)+il{t)) 



-Ve [x^{t) cos9{t) -x^{t) sin9{t);t] 

contains an explicitly time-dependent potential term 

, if z > and < R{t) < 00 
VE[z;t]^ E ^ 00 , if z < and < i?(t) < 00 

otherwise , 




(23) 



(24) 
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that vanishes at times t for which the classical path does not lie on the surface 



2.2 The Semiclassical Green function 

The semiclassical action of Eq. (|22(l can be interpreted as describing the non- 
relativistic motion of a particle of mass niE — E/v'^{E) in a homogeneous 
but time-dependent gravitational-like potential in two dimensions. To com- 
plicate matters there is a (generally time-dependent) restriction that Xr{t) = 
Xj.(t) cosd{t) +x^{t) sm9{t) > whenever < R{t) < oo. In analogy to motion 
in a gravitational field, the restriction Xr{t) > will be referred to as due to 
a (time-dependent) ceiling. Since the force is always directed away from the 
ceiling, there is no classical trajectory that starts from the ceiling and ends 
at the ceiling after a finite time t(7c,-E). However, the quantum mechanical 
amplitude for again observing the particle at the ceiling after a time t{jc,E) 
does not vanish. We now proceed to calculate this amplitude and thus describe 
diffraction seniiclassically. 

To simplify the notation, let us introduce the two-dimensional position- and 
momentum- vectors 

x = {Xj.,Xj_), ^= {mEX^,mEi±) (25) 

to describe the motion of the particle in phase space. The semiclassical Hamil- 
tonian H^^ of our two-dimensional problem depends on the classical path 7c 
and is explicitly time dependent. It is given by a Legendre transformation of 
the Lagrangian of Eq. (|23l) and in the notation of Eq. I|25|) reads 

^2 

(x,n;t) = K- k - = 7^ \- VE(xr(x,t);t) - E. (26) 

2mE 

Note that the parameter E of this two-dimensional Hamiltonian also determines 
the zero of the energy scale. This constant is unimportant for the dynamics at 
a given value of E, but is relevant when comparing the phases of the ampli- 
tudes at different wave numbers. Upon quantization of the classical problem 
described by the Hamiltonian of Eq. (|26|) , the amplitude that a particle in this 
two-dimensional world found initially at the transverse deviation x = a from 
the classical path will appear with transverse deviation x — b a. time r(7c, E) 
later is the causal Green function 



g^^{d,b;T{jc,E)) = (&|exp-^--^/ dtH^^{x,n;t) 



= y^[da;]exp^5-(7) . (27) 

The (time-ordered) exponential in Eq. (|27|l is the time-evolution operator gener- 
ated by a time-dependent Hamiltonian for a fixed time interval - here t(7c, E), 
for given E - and the paths 7 in the latter representation of the Green function 
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as a path integral in coordinate space are described by their transverse devi- 
ations, a at t = and b at t = T{'-fc,E), from the classical path 7c. (We use 
the notation Q and G for exact and semiclassical Green functions. Rather than 
introduce a new symbol, we choose to write Q^^ for the exact Green function to 
the semiclassical Hamiltonian of Eq. (|25|l for the classical path 7c.) 

The semiclassical energy Green function G'(X, Y; E) is proportional to the 
sum over all classical trajectories 7c of the amplitudes, J2-y^Sfc{^ — b = 
0; r(7c, E)). As noted in Appendix A, the two Green functions are related^) by 
a factor that depends only on the local group velocities VgCX-^E) and Vg(Y,E) 
at the endpoints of the classical trajectories, 

^' ^) ^ / L ,v E (« = = 0; r(7c, E) ) . (28) 

lhy^Vg{l^,E)Vg{Y,E) ^ 

For a homogeneous medium, the group velocity in fact does not depend on the 
location outside V and Wg(X, i?) = Vg{Y,E) — Vg{E) and Eq. (PSl) simplifies 
accordingly. 



3 An Example: Diffraction by a Sphere in a 
Uniform Medium 

To better understand the general case of Eq. (|27() , it is illustrative to consider 
the example of diffraction when the excluded volume V is a sphere of radius R. 
The exact expressions for the Green functions are known for this case and their 
asymptotic expansion for large kR = ER/{hv{E)) has been obtainedQl |S1 El- 
Besides making contact with previous work, the spherical case also exhibits 
some analytical properties of the Green function Q^^ that are essential in the 
construction of the solution for a more general obstacle. 



X 


X 




\ 


t=0 


\ X U '^x' V \ 


Y 



Fig. 2:The classical path, from X at t = to Y at f = Tc with winding number 
w = ^r\, for a spherical excluded volume V of radius R. The path is planar and 
and R fix the location of X with respect to the center of the sphere, 0, and lo and ly 
then determine Y. Note that for \w\ > 1, the classical path would completely encircle 
the sphere. The phase velocity is v{E), and Tx = Ix/v{E),To — lo/v{E),Ty = ly/v{E), 

with Tc = Tx +To +Ty. 
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The initial point X, the final point Y and the center of the sphere V 
define a plane. Every classical trajectory lies wholly in this plane and the angle 
a{t) vanishes. As is evident from Fig. 2, classical trajectories in this case can 
be assigned an integer w = 0, ±1, ±2 . . . that represents the number of times 
and the direction that the classical trajectory winds around the sphere in going 
from X to Y. w — Q corresponds to the direct trajectory and does not occur 
if, for given X, Y lies in the shadow cast by the sphere, w — ±1 are the two 
trajectories of opposite winding sense that do touch the sphere but do not wind 
fully around it, w = ±2 is assigned to the two trajectories that wind once but 
not twice around the sphere, etc. ... In addition to the sense in which it winds 
around the sphere, a classical trajectory 7c is characterized by three lengths, 
lx,lo and ly. Ix is the length of the straight line segment of the trajectory 
that is tangent to the sphere and extends from X to a point Xq on the surface 
of the sphere. Similarly ly is the length of the straight line segment of the 
trajectory that is tangent to the sphere and extends from Y to a point Yo on 
the surface of the sphere. Finally lo denotes the length of the classical trajectory 
that "creeps" along the surface of the sphere between Xq and Yq. The time 
intervals corresponding to these lengths are Tx = lx/v{E), Ty = ly/v{E) and 
To = lo/v{E), respectively, with 

Tc = t(7c, E) ^ Tx +To + Ty . (29) 

It is not difficult to ascertain that such a trajectory is of extremal length. 

The Green function G-y^ia = & = 0;Tc) is a function of the parameters 
E,v{E),R and of the variables Tx,Ty,To describing the classical path. Since 
a{t) = 0, it follows that — x,. and that x^ — Xi_. The potential in Eq. 
therefore does not depend on x±_ . The hamiltonian of Eq. H2t)|l in this case is of 
the form 

H^^^H^{TT^) + Hr{xr,'n:r;t)-E , (30) 

and the dynamics of the perpendicular and radial degrees of freedom separates. 
The time-independent Hamiltonian 

^ (31) 

describes the one-dimensional motion of a free non-relativistic particle of mass 
ruE = E/v'^{E), whereas 

HriXr,Trr;t)^ -^+VE{Xr;t) (32) 

ZiriE 

governs the radial dynamics. The boundary conditions a — b = imply that 
Xr{0) — Xr{Tc) — and x^{Q) = x^{tc) = 0. By virtue of the form of Eq. (|^ . 
the Green function Q-^^ {a — b — 0; Tc) decomposes into the product of two Green 
functions and an exponential, and we can write 

(0, 0; Tc) = exp{iETjh)go{0, 0; rc)e.(0, 0; r^) , (33) 
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where, for arbitrary x, y and r, 



^' = V » I — — ) ^ ^ 

is the Green function for a free non-relativistic particle of mass rn^ in one 
dimension and 

Gr{x,y;T) = {y\cxp 1^-^ J dtHr{x,7T;t)^ \x) (35) 

is the Green function corresponding to the time-evolution generated by the 
radial Hamiltonian of Eq. H32I) . The exponential factor in Eq. H33I) originates 
in the —E term in Eq. H3U|I and is the phase associated with the classical path. 
The explicit time dependence of the radial Hamiltonian in the case of a sphere 
is quite simple since the curvature of a classical path 7c is of the form 

^ _ f 1/R , <t <T^+To ^26) 



R[t) \ , otherwise 

In the time intervals [0, t^] and [tx +To, Tx + To + Ty] the Hamiltonian of Eq. H32|l 
describes the time evolution of a free non-relativistic particle of mass ttie- If 
the radial coordinate states were complete, Qr could be further decomposed as 

Gr{0,0;Tc) = da dbgo{0,a;Tx)gR{a,b;To)go{b,0;Ty) 

k f°°da db ^ , , , /"ika^ ikh' 
= Tr~ ~n= ^=&fl(a7o;T-o) exp 



^TJ"* JO Jo V 2^2; 2/j, 



(37) 

where we used Eq. H34(l in the latter expression. Here a and b denote the 
radial deviations of 7 from 7c at the points Xq and Yq of Fig. 2. We will see 
in Sec. 4.4 that the propagation to a point at the surface of the obstacle is 
only approximately described by the free Green function. (Roughly speaking, 
a photon with large but finite k effectively penetrates a small distance into the 
sphere.) The assumption that Qr separates in the form of Eq. (|37|l neglects 
subleading terms of the asymptotic expansion in the penumbra. 

g R is the coordinate representation of the time evolution operator generated 
by the time-independent Hamiltonian 

i7,(.,.) = ^-|. (38) 

in the half-space x > 0. The integrals over a and b in Eq. (|37|l extend only 
from zero to infinity because gf({a,b;T) vanishes whenever either a or 6 are 
negative. We recognize that the one-dimensional quantum mechanical problem 
described by Eq. (|38|l is the propagation of a non-relativistic particle of mass 
iriE = E/v'^{E) under the influence of a constant force leading to an acceleration 
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\g\ = v'^{E)/R directed away from x — 0. Note that Qii{x,y;T) — ^^(y, a;; — r), 
because the Hamihonian Eq. (|38() . although not bounded below, is hermitian. 
By definition, the causal Green function Qn satisfies the second order differential 
equation 

f fi^ d'^ Ex\ d 

in the half-space x > 0, y > and the initial condition 

gRix,y;0) = S{x~y) . (40) 

The solution to Eq. (|39|l and Eq. (|40|) is completely specified by imposing the 
boundary conditions 

gR{0,y;T)^0 and lim ^^(x, r) = , (41) 

X — >OG 

and requiring that the magnitude of Qn remain bounded for t — > oo 

IGnix, y;T ^ oo)\ <oo . (42) 

The latter condition excludes any solution which increases exponentially and 
is just the requirement that the presence of an obstacle should not lead to the 
production of photons in the shadowed region. 



3.1 The spectral representation of ^ 

The Hamiltonian (|38|1 is hermitian and bounded below only for negative values 
of the parameter R. In this case it can be interpreted as describing a non- 
relativistic particle in a homogenous gravitational field in the half-space a: > 
above a table. The spectral representation of the solution to Eqs. H39I) and H40|l 
that satisfies the boundary conditions in Eqs. (|41|l and (|42(l is well known iff)). 
Our strategy is as follows. For R > waves diffract and there are no classical 
orbits, while for R < there are classical orbits, a much simpler situation. We 
will therefore analyze the R > case by first considering the R < case, and 
then return to the R > case of interest by analytic continuation. 
With R replaced by —R, we have 

k ^ 

Q-nix, y- To) = ^ 2_. ^n{x)^n(y) exp {-'iTov{E)e,,a) , (43) 

n=0 



where the inverse length scale a is 

k 



2l/3(yfci?)2/3 ' 



(44) 



and the wave function 



vI/„(.x)=Ai(^-6„)/|Ai'(-6„)| (45) 
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is given by the Airy function[TF| and its derivative, Ai'(z) — {d/dz)A\{z). Con- 
sidered as a function on the complex plane, the Airy function Ai(z) is the 
solution to the differential equation 



dz 



^Ai(z) ^ zki{z) (46) 



with the asymptotic behavior for \z\ —f oo, 

1 2 
Ai(z) - -^z-i/'*exp(--z'''/2) , |argz| < TT 

Ai(-z) ^ -L,-i/4 3in(2,3/2 ^ i^^g^i < I ^ (47) 

The first of the boundary conditions of Eq. H41|) relates the e„'s to the zero's of 
the Airy function 

Ai(-e„) = , 71 = 0, 1, . . . (48) 

Since the zero's of the Airy function are all on the negative real axis|18|. the e„ 
are real and positive and are ordered as < e„ < Cn+i- That Ai'(— e„) gives 
the proper normalization of the wave functions follows from the identity, 

dxAi^{x) = (Ai'(z))^ - zAi^{z) . (49) 

Eq. I|49(l is proven by taking the derivative of both sides and using Eq. (|46|l . 
Choosing z — — e„ in Eq. I|49|l . dividing by [Ai'(— e„)]^, and redefining the 
integration variable as a; ax— en, one finds that dx'^'^{x) = 1/a. Since the 
are eigenfunctions of a hermitian operator corresponding to the eigenvalue 
e„, ^„ and are orthogonal for n =/= m and therefore. 



f 

Jo 



dx'^nix)'^m{x) = Snm/CT . (50) 



The completeness of the spectrum of the bounded hermitian operator H_ji in 
the space of normalizable functions on the half space 2: > that vanish at x = 0, 
together with Eq. proves that Eq. igSl also fulfills Eq. Thus the 

spectral representation Eq. (|43|l would be the (unique) solution of Eq. H39|l 
and Eq. H40|l we are looking for if the parameter R in Eq. H39|) were negative. 



3.2 The spectral representation of Qr by analytic contin- 
uation 

To obtain the spectral representation of Qj^ for positive real values of R, we 
analytically continue Q-r in the complex plane. With R the modulus of the 
complex number i?^, 

i?0 = Re"^ , (51) 



15 



the analytic continuation of the spectral representation Eq. (|43|l G^r^ in the 
phase (j) is uniquely given by requiring that Eq. H42() hold for all < < tt. 
The dependence of the inverse scale factor a on the phase is 



ere 



(52) 



For sufficiently large positive values ofx, < a,rg{kx/{a^R^) — en) < for any 
fixed n and the Airy function decays exponentially as given in Eq. 1)47(1 . Note, 
however, that the behavior of the wave functions for large arguments becomes 
oscillatory with a power-law decay at the endpoint (j) = ir oi the interval. This 
is consistent with the fact that there is only one turning point of a trajectory 
with fixed energy in this case. The boundary conditions Eq. ((41|l thus hold for 
all (j) £ [0,7r]. Finally, because 

Im(T0 < for < < TT , (53) 

the Green function of Eq. 143|) decays for large times r and Eq. 1)42(1 holds. Note 
that the only dependence of Qji on the arc length Ig of the creeping section 
of the trajectory is via r = Tq = lo/v{E). thus decays exponentially with 
the length Iq of the creeping section with an exponent that is proportional 
to loikRY/^ / R. This dependence of the exponent on the cube root of the 
wave number is somewhat unusual and perhaps unexpected for a semiclassical 
expansion, but is typical for diffraction^. Setting (/> = tt in Eq. 1(43(1 . R R^, = 
—R and therefore cr — > (Ttt = ct, where 



2l/3(fci?)2/3 ■ 

Using Eqs. ((43(1 and 145(1 . one thus finally obtains 



(54) 



(55) 

as the spectral representation of the Green function that satisfies Eqs. ((39(1 
and ((4U(I as well as the boundary conditions Eqs. ((41(1 and ((42(1 . 

Inserting Eq. ((5^ in Eq. (|S7I) . Gr is found to have the representation 

Gr{0, 0; Te) = J2 ^" ('^^-) ^» (^^2/) ' (56) 



AnaR 

n=0 



with amplitudes 



^ r, Ai(pV?-e„) fip^\ 

""-^^^^L |Ai'(-e„)| '^^H-J ' ^''^ 

the result previously obtainedj^j for the asymptotic behavior in the deep shadow 
region. Within the semiclassical approximation it might appear to be consis- 
tent to further approximate the Airy function and its zeros semiclassically and 
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evaluate the integral in Eq. H57|l in saddle point approximation. This is done 
in Appendix B and indeed was part of the original approach!] 21- However, al- 
though perfectly reasonable, one should perhaps point out that it is not simple 
to justify this procedure from the semiclassical point of view: the error of the 
semiclassical approximation is largest for the lowest zeros of the Airy function, 
which apparently dominate the sum in Eq. H5t)|l in the semiclassical limit. An 
alternative to the semiclassical approximation of the spectral representation of 
the Green function and associated integrals is the direct semi-classical evalu- 
ation of the path-integral discussed below. It gives corrections to the above 
procedure that become important outside the deeply shadowed region. 

4 Diffraction by a Sphere: Semiclassical Evalu- 
ation of the Path Integral 

In the previous section we constructed the exact spectral representation, given 
by Eqs. (|55|) and (|57|l . of the Green- function Eq. Ij^Sfl for the radial motion 
by analytic continuation. At the end of the calculation we argued that the 
semiclassical approximation of the zeros, Cn, and of the remaining integrals I?„ 
may not be fully justified. We now proceed to do the calculation in reverse; 
we first obtain the semiclassical approximation to Gr for negative values of the 
parameter R and then analytically continue the semiclassical result to positive 
values of R. Although this procedure is of only pedagogical value in the case of 
a sphere it has the great advantage of being generalizable to more complex 
situations. We will see that this method of evaluation is also interesting for 
conceptual reasons. 

The point is that there are classical trajectories for the radial motion from 
a; to 2/ in time r if the parameter R in the definition Eq. (|24|) of the potential 
VE{Xr',t) of the Hamiltonian Hr given in Eq. (|32|l is negative. In this case the 
force on the particle is either absent or directed toward the endpoint a,t Xr = 0. 
The force of magnitude E/R\s turned on at "time" Tx and lasts only for a finite 
time To- Using the horizontal axis to denote the time and the vertical one to 
denote the "height" x{t), the classical trajectory during the time interval To is 
that of a ball bouncing up and down vertically on a table. See Fig. 3. There are 
many classical trajectories that start with an initial height a;(0) = x and which, 
at Tc = T^j-t-To-f Ty, attain a final height x{tc) = y. They are distinguished by the 
number of bounces "SO", n = 0,l,2,.... If the initial or final point are above the 
table, the classification of a classical trajectory by the number of bounces is not 
unique. There are classical trajectories with the same number of bounces that 
differ in the directions of the initial and/or final momentum. The ambiguity for 
a given n is at most 4- fold. If a; < and y < 0, this ambiguity is removed since 
the ball has to be above the table in the intermediate region. The signs of the 
initial and final momenta must be positive and negative, respectively, and the 
classification of a classical trajectory by the number n is unique. Eventually we 
are interested only in trajectories that start and end level with the table. To 
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avoid unnecessary (and irrelevant) complications, we consider only the situation 
where the initial and final points satisfy x < and y < 0, that is, are below or 
level with the table. 



x.(t) 













n=2 







Fig. 3: Radial deviations with n = (no bounces) and n — 2 (two bounces) from a 
particular classical path 7c for a spherical obstacle of radius R. In the time intervals 
[0, Ta;] and [tx + ro,Tc], there is no acceleration and these sections of the paths are 
straight lines; the motion in the interval [tx,Tx + To] is governed by a constant accel- 
eration g = —v'^{E)/R. The initial and final radial deviations from ■jc, x — a;„(0) and 
y = Xn{Tc), not labeled in the figure, do not depend on n. t^"^ is the i-th turning point 
of the path with n bounces, r^"' is the time between successive bounces, or successive 
turning points. 

The classical trajectory 7„ connecting x and y in time with n bounces is 
a piecewise connected solution to the equation of motion 



Xnit) = -.g(t) 



(58) 



with x{0) — X and x{tc) — y and x{i) > for < t < + Tq that apart 
from the n points corresponding to bounces is continuous in phase space. Since 
the (negative) acceleration g{t) = v'^{E)/R — g does not depend on time in the 
interval r^, < i < Xj, + Tq, the "energy" (for notational convenience we extract 
the scale to have a dimensionless energy) 

is a real positive constant of motion that characterizes the trajectory in this time 
interval. In the initial and final time intervals t € [0, r^,] and i € [t^, + To,t^ 
the acceleration vanishes and the momenta are constant. We also introduce 
dimensionless initial and final "momenta" , 

Pn = TTn{0)/{nk) = 7r„ (t^, ) /(?ifc) > , 

Pn = -TTn{Tc)/{nk) = -7r„(r, + T^) / {hk) > . (60) 

The signs have been chosen so that p„ and p'j are both positive for initial and 
final points x < and y < 0. The equation of motion l|58|) with g{t) = implies 
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that 



Xn{Tx) ^ X + IxPn and Xn{Tx + To) = y + lyp'n . (61) 

Demanding continuity of the classical trajectory in phase space at t = r^; and 
t — Tx + To, we determine p„ and p'^ in terms of e„. They are given by 

Pn = Ve„ + [Ix/RY - 2x/R - Ix/R 

p'^ = ^en + {ly/RY-2y/R~ly/R. (62) 

The signs of the square roots have been chosen so that both scaled momenta, 
Pn and p^, are positive for initial and final heights a; < and y < 0. 

To completely specify the trajectory, we need an equation for e„. The ac- 
celeration is constant for < t < + Tq and the first turning point of the 



trajectory occurs aX t = r}"'' , where 



t; = Tx 



-PnR/v{E) . (63) 



The last turning point, the (n + 1)"^, occurs at t = where 

Ti%=Tx+To~p'nR/v{E) . (64) 

The equation of motion (|58|l for constant acceleration finally relates the time 
r^"^ between successive bounces (or successive turning points) when there are 
n bounces in all to the conserved "energy" e„, 

4"^ = 2^R/v{E) . (65) 

The time to the first turning point plus the time for n bounces plus the time 
from the last turning point to the end of the trajectory is the total time Tc of 
the trajectory. We thus have that t}"-* + nr^g'' + (jc — Tjj'J_\) — Tc, or 

Pn + Pn + ^Uy/E^ ^ lo/R . (66) 

Eq. (|66|l together with Eq. (|62|) is a fourth order algebraic equation for e„ whose 
roots can be found analytically. There is only one real and positive solution e„, 
because the left hand side of Eq. (|56|l is a monotonically increasing function of 
e„ that vanishes at e„ = and goes as 2{n + for large values of e„. We 

remark that en,Pn and p'^ given by Eq. H66|l and Eq. (j62|l depend only on the 
lengths lx,lo and ly of the problem and the scaled initial and final heights x 
and y (all measured in units of R) . The solution e„ of Eq. (|66|l therefore does 
not depend on E nor v{E). The action Sn{x,y) of the classical trajectory with 
n bounces is proportional to ER/v{E) — TikR. After some algebra, Sn{x,y) 
can be expressed in terms of CmPn and implicitly given by Eq. H66|l and 
Eq. (inSl. We find 

Sn{x,y) ^ Pn + Pn + 2nef/^ p^Jx +Pnly ^ sjo 

hkR 3 2R ' ^ ' 
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Using Eq. H62I) . one explicitly verifies that 



dSn(x,y) hkR , . , , , ^ , , 

—ff^ = —(Pn+ Pn + - lo/R) . (68) 

Eq. thus selects the energy e„ at which the action is stationary. This fact 
will be of some consequence when we compare the present approach to the one 
of section 3. 

The semiclassical approximation Gr{x,y;Tc,~R) to the Green function 
Grix,y;Tc, —R) defined in Eq. H35|l might seem to be the sum of the asymp- 
totic contributions from the infinite number of classical path. The semiclassi- 
cal contribution of a single classical path is known; it is given by Van Vleck's 
formula|21j. If we could simply sum over the asymptotic contributions of all 
classical trajectories, we would (in this one dimensional case) obtain 



exp I 

x—y—O 



Sn{0,0) 37m\ 
h 2~J 



(69) 

(There are n conjugate points and n points of reflection on the classical trajec- 
tory with n bounces, contributing in total a phase factor of exp(— i37rn/2).) 

Although each term of the sum is the correct asymptotic contribution from 
a particular classical trajectory, we now show that Eq. H69|) does not give the 
asymptotic behavior of the Green function Q,. for kR oo. The cause of the 
problem is that for any fixed value of kR, the classical trajectories and their 
corresponding action S'„ become indistinguishable as n —^ oo. The "floor" 
at a; = is a zero-dimensional caustic of the classical motion^^. It then is 
not possible to interchange the implicit kR —^ oo limit with the sum over n 
in Eq. P^ . 

To see that the interchange is not allowed, we examine the asymptotic be- 
havior of the classical action S'„ for trajectories with large n, assuming for the 
moment that Eq. H69|l is valid. We begin by observing that the second derivatives 
of the action appearing in Eq. H69(l can be evaluated in terms of the solutions 
Pn,Pn and e„ of Eq. and Eq. for x — y = 0. The classical equations 
of motion imply that the initial momentum 7r„ (0) is given by 

^^H^ = -7r„(0) = -nkp^. (70) 

Eq. (|70|) can be checked explicitly by taking the derivative of Eq. 16 7|) and 
using Eqs. and (|^ . From Eq. T7U\i we thus have that 



d'^Sn{x,y) 



dxdy 



= nk'-^ 

x=y=o dy 



(71) 

x— y— 



The variation of the initial momentum with the endpoint of the trajectory re- 
quired in the latter expression is implicitly given by Eq. H66f) and the defini- 
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tions Eq. l()2|l . After some algebra we obtain 



(72) 



Note that in Eq. (|69|l we need only know the right hand side of Eq. (|72|l at 
X — y = 0. To evaluate the function Eq. (|69|l it thus is sufficient to solve Eq. H66|l 
with X — y ~ 0. 



The solution e„ of Eq. (|66|l will be small compared to either l^/B? or ly/R^ 
when n is sufficiently large. For x — y — 0, we then may expand the roots 
occurring in Eq. Ht)2|l and thus solve Eq. (|66|l asymptotically for large n with 
the result 

= ^ (lo/R) + 0{n~\x/l^,y/ly) . (73) 
The scaled momenta p„ and p'^ defined by Eq. Ht)2|) therefore are of order n 



-2 



and 

d^Sn{x,y) 



dxdy 

Note that in this limit 



x—y—Q 



flklo 



2 



0(n-4) . (74) 



depends neither on nor on ly. This is readily understood. For n bounces 
within a given time interval Tc, their height becomes exceedingly small compared 
to Ix or ly. For x = y = and large n, the initial and final momenta of the 
trajectory are of order and the free sections of the trajectory therefore do 
not contribute appreciably to the classical action. One can check that Eq. H75|l 
indeed is just the action for a classical trajectory with n bounces that altogether 
take a time Tq = lo/v{E). The asymptotic behavior in Eqs. (|74|l and H75() implies 
that the sum in Eq. I(j9|l converges to a logarithm for large n. 

The asymptotic form Eq. (|75|) . however, also shows that the actions of clas- 
sical trajectories differ by much less than h for n ^ y/kR. For a fixed value of 
kR, the classical paths with sufficiently large n therefore could (and should) be 
considered quantum fluctuations of each other. It becomes a matter of exchang- 
ing limits: the asymptotic weight of a particular classical path characterized by 
n bounces in the limit kR — > oo is indeed given by its contribution to Eq. (|69|l . 
However, we actually would like to know how the (infinite) sum of classical paths 
contribute to the Green function for a large, but fixed value of kR. The problem 
arises because the trajectory lim„_>oo a;„(t) = a;oo(^) = is a caustic that is 
approached by an infinite number of classical trajectories. Note, however, that 
the sum in Eq. H()9|l would be cut off at a finite value n = rimax for a negative 
value of the initial point x or final point y. Due to the non-vanishing initial 
and/or final momentum, the "energy" e„ of a trajectory in this case is bounded 
from below, e„ > Max ((x/^^r)^, (y/ly)^), and Eq. can be solved only for 
n < f^max- Eq- H69() in this case would give the correct asymptotic expression 
for sufficiently large kR. 
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For a given geometry and wave number k, successive terms of the sum 
in Eq. (|69|l tend to interfere destructively for sufficiently high n - one thus 
should reorder the sum and collect terms of similar phase. To do so, we make 
use of our earlier observation in Eq. H68|l that Eq. H66|l gives the stationary point 
of Sn with respect to a variation of the "energy" e„. The second derivative of 
the action at this stationary point is 



hkR 



del 



Pn+Pn- lo/R 



Pn + Ix/R P'n + ly/R 



hk[loiPn 



l){Pn + '-^)+PnUPn + ^) + P'JyiPn + 'f)] 



Aen{pn+lx/R){p'n+ly/R) 



x—y—0 

> 



x—y—0 



(76) 



where we have used the definitions (|62|l to write the numerator in terms of p„ 
and p'„ only. 

We will show just below that it is consistent, within the semiclassical ap- 
proximation, to replace Eq. (|69|l by the integral expression 



a(0,0; re,-i?) 



where 



e de B{e)e 



iS(e)/n 



n=0 



,inf{e) 



2kR 



(77) 



(78) 



In Eq. fTfl) the factor B{e) is given by Eqs. ifT^ and lf75|) with e„ replaced by 
e^, that is. 



B{e) 



dxdy 



del 



hk 



x—y—Q;en — e'^ 



2[(e2 + (/./i?)2)(e2 + (/,/i?)2)]V4 

and ^(e) is related to the action of Eq. (|67l) . 

p^l^+p'% ~ e% 



5(e) = hkR 



2R 



(79) 



(80) 



where here too e„ has been replaced by e^. In addition the only explicitly n- 
dependent term in Sn{x,y) of Eq. has been separated out. In Eqs. ifT^ 
and (|80|l the variables p and p' are the dimensionless initial and final momenta 
of Eq. for X = y = and "energy" e^, that is 



P=V^^TJlJW-lx/R, P' ^^e^ + {lylRf~ly/R 



(81) 



To show the semiclassical equivalence of Eqs. (|69|) and (|77|) . recall that the 
phase of the integrand in Eq. (|77|l for a given n is stationary at e = Je^- 



22 



Expanding about e = y^e^, and evaluating the resultant Gaussian integral, one 
finds that Eqs. H69|l and (|77|l are indeed equivalent term by term. 

The asymptotic contribution of a classical path with a given n to the Green 
function in Eq. H77|) thus is precisely the same as in Eq. (|69|l . On the other hand 
we now can interchange the limits and obtain the correct asymptotic expansion 
by summing the geometrical series in Eq. (|77|l for a given value of kR. With a 
small positive imaginary component of e, the sum in Eq. (|77|l is evaluated as 

oo 

lim e™(-^(^)+"') = lim tjtt^- (82) 

n—O 

The integrand of Eq. (|77() therefore has simple poles just below the real axis 
at /(e) = 27rm — ir] for integer m. Apart from these simple poles, the integrand 
is analytic in the shaded region of Fig. 4. The integration over e in Eq. H77|l 
along the real axis therefore is equivalent to integrating along the contours Csd 
and Cji shown in Fig. 4 and including the residues from the enclosed poles. The 
paths Csd and Cn in the complex plane are chosen so that, as shown below, 

i) the contribution from C-jz to the integral is negligible for sufficiently large 
TZ, and 

ii) the modulus of the integrand in Eq. H77() decreases monotonically and as 
fast as possible along Csd, the path of steepest descent. 

The semi-classical evaluation of the integral in Eq. (|77|l thus has three distinct 
contributions, to be denoted by Gpoics, Gctz and Gc^^. We discuss them sepa- 
rately. 



Im e 




Fig. 4: Contours of integration. Apart from discrete poles just below the positive real 
axis, the integrand of Eq. I77II is analytic in the shaded region of the complex plane. 
The integral along the positive real axis in Eq. 1771 from the origin to TZ is evaluated by 
integrating along the contours Csd and Cn, taking the residues of the poles inside the 
closed contour into account. Csd is the path of steepest descent of the integrand that 
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starts at the origin. It asymptotically approaches the dashed line with arg e — — 7r/4 
for large |e| oo. 



4.1 The Contribution from Poles with Integer m > 

In the shaded region of the complex plane of Fig. 4, the integrand of Eq. H77|l 
is analytic except for simple poles at 

3^(m + 3/4)^ ^ ^ ^ ^'^ 



just below the real axis for any integer m > 0. Note that Cm corresponds to a 
semiclassical determination of the (m + l)**^ zero e™ of the Airy functionjl6|. 
The contribution Gpoios(— i?) to Gr(0, 0; Tc, —R) from the poles is 



oo . 



residue of \ ^ ^ B{e{m)) ^^s(e(m))/h (-g^^ 



pole at e(m) J ^ifci? ^ e(m) 

m— ^ m— 

Using Eqs. (O lO, (jHIJ), (EOIl and Gpoics(--R) can be rewritten as 

k 

n=0 

with diffraction amplitudes 

The Dn{£,) depend on e„ defined by Eq. and the corresponding (rescaled) 

" momenta" 

Pn{S.)^^en+e-^ ■ (87) 

Since p„(^) ^ e(n) ^ n^/^ for large n, the n*'' and [n + 1)*'^ contribution 
to the sum in Eq. H85|) differ by a /inzie phase. By first summing in Eq. H77|l 
rather than semiclassically evaluating the integral over e for each summand, one 
avoids the problem discussed earlier that for fixed wave number fc, classical paths 
near a caustic interfere with one another. [Note that the individual summands 
of Eq. H85|l do not correspond to contributions from individual classical paths 
of the original problem.] 

The similarity between the expressions in Eqs. (|85|l and (|56|l is no coin- 
cidence. We show in Appendix B that apart from approximating the zero's 
of the Airy function e„ by e„, is proportional to the integral of Eq. H57|l 
evaluated in saddle point approximation. We thus could have obtained the ex- 
pression Eq. I|85|) by evaluating all integrals in saddle-point approximation and 
replacing the wave- functions 5'„(a;) and corresponding energies e„ in Eq. H56|l 
by their WKB counterparts. The fact that such a procedure ignores any asymp- 
totic contribution from the integration along Csd and C-n indicates that this may 
not be the whole story, and indeed it is not. 
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4.2 The Contribution 

The contribution to the Green function from integrating along the arc C-ji is 
neghgible for sufficiently large TZ. To see this, we use Eqs. ((77|) - to write 
this contribution to Gr in the form 

i r i[{S{e)/K)-f{e)\ 

GcA-R) = ^ J^^edeB{e)^—-j^ , (88) 

where the contour C-n is shown in Fig. 4. 

Consider the behavior of the integrand in Eq. H88|) . with the contour C-n 
chosen to be an arc of very large radius TZ — |e| centered on the origin that 
extends from the positive real axis into the lower half of the complex plane with 
— 7r/3 < arg e < 0. The real part of if{e) on the arc is positive and proportional 
to |ep. e^^^^"^^ therefore is exponentially small on Cu and for sufficiently large 
TZ = |e| is negligible compared to 1 - the denominator of the integrand thus 
approaches unity. To expand {S [e) / fi) — f [e) for large |e| and — 7r/3 < arg e < 0, 
note that p '-^ e — Ix and p' ^ e — ly vo- this region. The leading term of 5(e) 
for |e| oo is cancelled by the leading term of /(e) and one has, for |e| ~ oo 
and — f < arg e < 0, 

{S{e)/h)~f{e)^~^{lx + ly + lo)e^ . (89) 

Since Re(— ze^) < 0, Eq. H89(l shows that the contribution of C-n to the integral 
in Eq. 188|) vanishes in the limit TZ ^ oo. 

4.3 The Contribution Gc^^ 

The contribution Gc^^{—R) is given by Eq. H88|) . with C-n replaced by Csd- The 
section Csd of the contour is chosen to coincide with the path of steepest de- 
scent that begins at the origin. From Eq. we see that arg e along this path 
approaches — 7r/4 for very large values of |e|. Because p and p' are both propor- 
|e|, the leading behavior in this limit is (using Eqs. H78|l 



and (IHOJ)) 



(5(e)A)-/(e) ^ - 



2 

-| < arg e < . (90) 



— -e -I e 



Oie% 



When lo > 0, the path Csd of steepest descent starts at the origin of the 
complex plane with arg e — — 7r/4. Unless R is very large compared to Ix 
and ly, the terms of order e^ of Eq. H9U|) are negligible in the semiclassical 
approximation. However, at grazing angles, lo/ R ~ 0, the leading dependence 
on the wave number comes from the term of order in Eq. (|90l) . To obtain a 
semiclassical expansion that is uniformly valid in the whole classically shadowed 
region, one therefore has to retain all exponential terms with exponents that 
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are of cubic or lower order for small e. To semiclassical accuracy we have, 
using Eq. a given by Eq. EH), B{0) by Eq. ^ and /(e) by Eq. 

Gc.A-R) = -T / edeB{e) — 



e 



27rZ?i(fci?)2/3 \^\^4y °y 4TTiy^ \\4: 

(Note that this contribution to the Green function neatly separates into the 
product of free one-dimensional Green functions over distances Ix and ly and 
a factor that depends on Iq and R.) For positive real values of z, the function 
P{z) has the integral representation, 

P{z) = im Cd£ . (92) 

P{z) is an entire function that is defined in the whole complex plane by the 
convergent series (derived in Appendix C), 

n = ) - (I) (---)•'" ^-.■./3(-.), W 

where 

°° n 

Li^(p) - E (94) 

is the poly-logarithm of order u. Of particular interest are the asymptotic values 
of P{z) for large and small z, corresponding to large and small values of klo'- 

1/3 

0.966e"°"347- 



^(0) - r(2/3)Li2/3(-z) 



~ ~^776 ~ ; , for |2:| - c», -- < arg z < - . 

(95) 

The asymptotic forms of Eq. (|95|l are most easily obtained from Eq. (|92|l : see 
Appendix C. 

To finally obtain Gr(0, 0; r^, i?) we again analytically continue the expres- 
sions in Eqs. (|91|l and (|85|l to negative values of R in such a manner that Gr 
remains bounded for large values of lo- We thus find that the semiclassical 
approximation to the Green function for the radial coordinate from a single 
extremal path whose end-point is in the shadow of a sphere of radius R is 



G,(0,0;T„i?) = *^ 



P ((3/4)1 a/ 
22/3a./IX 



(96) 

The diffraction amplitudes £'n(C) a-rc defined by Eqs. (|86() and H87I) . with e„ 
given by Eq. I|83() and a by Eq. (|54() . The function P{z) for complex z is defined 
by the expansion Eq. (|93|1 . 
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4.4 The Penumbra near the Glancing Ray with lo = 

The correction from the contour Cgd in Eq. (|96|1 is essential for describing the 
penumbra. We compare our result with the asymptotic form of the exact solu- 
tion in the penumbra obtained injSlin|. In the penumbra klg <C kR and one can 

expand Gr in powers of klo- The leading term in this expansion of Gr is 



G.(0,0;re,i?)|, 



ik 



AnaR 



P{0) 



Dn {crlx) Dn (crly) 



V n=0 



(97) 



The sum in Eq. H97I) is exponentially damped and converges for any value of 
kR. Let ^ — al, with / — Ij. or ly\ a is given by Eq. H54|) . For e„ 3> we find 



on using Eq. that 
and therefore that 



(98) 
(99) 



Since Im(^) < 0, -Dn(C) defined by Eq. H8t)|l therefore decays exponentially and 
the contribution to the sum in Eq. H97|) from terms with e„ ^ is negligible. 
We now assume that e„ « For |e„/^^| « 1, 



(100) 



Since |^| oc (fci?)^/'^, we can simplify the sum in Eq. H97() in the asymptotic limit 
we are interested in. Retaining only leading terms in kR, the definition Eq. H86|l 
of the diffraction amplitudes leads to the simplification 



E_D = ^ Dn {(Tlx) Dn {aly) 



n=0 



exp 



. t„ la; -r ly 



4ct Ixly 



LXJ 

^^$(n) (101) 



The sum over n in Eq. HlOlfl is effectively cut off when becomes of order 1/ A 
or larger, where 

_ f-x + ly 



\a\lxL 



(102) 



Since for large values of ^, e„ ^ 1/Lanibda ~ |^| « the sum in Eq. (|101|l 
effectively never extends to values of n where the expansion of Eq. H1UU|) is not 
justified. 

We appeal to the Abel-Plana formula P^. 



n=0 



dn $(ri) + i 



dz 



$(iz) - $(-iz) 
exp(27r2:) — 1 ' 



(103) 
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to evaluate the sum in Eq. puill asymptotically. In the evaluation of the first 
and third terms of Eq. H103() , we approximate the exponential factor in Eq. (|101|l 
by unity. We find, using Eq. that, 

^ ' - ' -I ' ^ + 0(A) , (104) 



2 2y% \9tt 

to leading order in A. To the same accuracy, the third term of Eq. H1U3|I becomes, 

Im[(iz+ |)-*1 

,^ t , + 0(A) . (105) 



3'''' / Jo exp(27rz) — 1 

To evaluate the integral over n of the second term in Eq. Hl()3|l . we change the 
integration variable from n to In, using Eq. (|83|l . We then find that 



dn^n) = /^^^^^^ * exp 



.6 Ix ~^ ^ii 
I — 

Act Ix^y 



ialjy (97r)2/3 



7r(/^ + ly) Air ' '^''^^ ■ 
(106) 

Collecting the results and using the them in Eq. (|97|l . the Green function at the 
grazing angle is found to have the asymptotic form. 



(107) 

For large values of kR, the Green function at the grazing angle is just half 
the direct term, with an additional contribution jSi, proportional to {kR)'s . The 
proportionality constant K of the latter is 



(97r/4)i/3 y i J \3 J Jo exp(27rz) - 1 

= P(0) - 2ReP(0) = -P*(0) « -0.966e~° ™3'^^''* . (108) 

{—)IleK differs from the coefficient c « 0.996 for the glancing contribution to the 
Green function of Rubinow and Wu [S] by about 3%; we find ImP(O) = —ImK = 
to be very small, while in|S] it is found to be zero. The small discrepancy arises 
because the integral over a ratio of Airy functions that defines c is evaluated 
numerically in rather than in saddle point approximation. Our approach is 
perhaps slightly more consistent - that does not imply that it is more accurate 
- and in any event, the difference is quite small. 

It perhaps is of some interest that the correction of order (fcP)^/^ for the 
glancing ray does not appear in Gr{0, 0; Tc). Eqs. and (pTzjl imply that 

e.(0,0;r,)L = -^1^1 ^<5(a-6)exp^*^"' 




(109) 
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without further corrections. Physically, this discrepancy can be traced to the 
assumed decomposition of the radial Green function in Eq. H37|l . The lower 
bound of the integral over the radial coordinates in Eq. (|37() can be set to zero 
only for fc = oo. For large but finite wave number fc, the correction of order 
(kR)^/'^ is reproduced by assuming that the fluctuations effectively penetrate 
the obstacle a small distance do{k, R). Replacing the lower bound of the radial 
integrals in Eq. (|109|l by —do{k,R), one finds, 





9 \/ 9 Y/ -u/ ^ o ■ /r^ + Oidt) . (110) 
Comparison of Eq. (|110|l with Eq. H107() leads to 

d„(i?,fc)^-jr ^^^^^^" (111) 

for do{k,R). The effective depth of penetration thus vanishes rather slowly, as 
^-2/3^ The phase e^*/^ in the expression Eq. (|111() for do has its origin in the 
analytic continuation of R to — i?, and the proportionality constant K is very 
close to —1. 

The true semiclassical Green function does not in general separate into 
free- and creeping- Green functions, and approximately reproduces the sublead- 
ing asymptotic contributions as well (to within about 3%). 



5 Discussion, Generalization and Conclusion 

The semiclassical description of diffraction is closely associated with extremal 
classical paths that satisfy Fermat's principle but are not stationary. Such paths 
arise due to the non-holonomic constraints imposed by "obstacles". The asso- 
ciated Lagrangian is given in Eq. (|23|l . It depends on the transverse deviations 
Xj.{t),Xj_(t) from the classical path and describes the two-dimensional motion 
of a non-relativistic particle of mass tue = E/v^{E) under the influence of a 
spatially constant but in general time-dependent, force mEg(t). The acceler- 
ation g{t) is inversely proportional to the curvature of the classical path and 
vanishes for time intervals in which the classical path 7^ is not constrained by 
an obstacle. On "creeping" sections of the classical trajectory, fluctuations are 
in addition restricted to the half-space exterior to the obstacle. 

Depending on the sign of the curvature, g{t) is either directed toward or away 
from the obstacle's surface - corresponding to the presence of (a generally time 
dependent) "floor" or "ceiling" in an analogous 1-dimensional gravitational- 
like problem. Motion in the presence of a "floor" is relevant for "whispering- 
galleries", for example, whereas the presence of a "ceiling" corresponds to 
diffraction into the classical shadow of an obstacle. 
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Diffraction is semiclassically described by tlie ampiitude tlrat the particle 
moves from tlie ceiling to the ceiling in a fixed time interval with acceleration 
g{t). This motion is classically forbidden. To calculate the amplitude semiclas- 
sically, we relate it to the one for motion from the "floor" to the "floor" with 
acceleration g{t). The latter problem, that of a "whispering gallery", admits 
classical solutions. The amplitudes of the two problems are related by analytic 
continuation in the phase (f> of the curvature R{t) |_R(t)|e*'^ of the classical 
path. 

The existence of a caustic, the "floor", that is approached by an infinite 
number of classical paths poses an additional problem in the semiclassical de- 
scription. For any fixed value of k, there are distinct classical paths whose ac- 
tions differ by much less than h. The fluctuations about different classical paths 
therefore cannot be separated and tend to interfere with each other. Attempt- 
ing to evaluate the path integral asymptotically by summing the asymptotic 
contributions of every individual classical path becomes highly inaccurate. 

The problem was solved for the special case of a classical ray that partially 
creeps along the surface of a sphere; the asymptotic analysis of the same problem 
was also at the heart of the original geometrical interpretation of diffraction ^E]. 
The curvature of this trajectory is piecewise constant and corresponds to an 
acceleration g{t) that is piecewise constant in time. A path with n bounces in 
total time Tc is stationary with respect to a variation of the (conserved) "energy" , 
e„, of the particle in the region of constant acceleration. One thus can express 
the asymptotic contribution from a classical path to given n as the result of an 
integral over the "energy" e that is evaluated semiclassically at the saddle point 
e — Bn. This manipulation is consistent with the semiclassical approximation 
and allows one to sum over n for a given "energy" e before evaluating the 
remaining integral over e asymptotically. 

In the case of a sphere, the exact expression for the Green function of a mass- 
less particle is known and has been evaluated asymptotically in the umbra^ of 
the sphere as well as in its penumbra|niini- Quite different asymptotic forms were 
found for the two cases. The procedure outlined above and described in more 
detail in section 3 gives a uniform asymptotic expansion that is valid in both 
regions: the exponentially decaying terms arise from pole contributions to the 
final integral over e, whereas power corrections (previously obtained separately 
for the penumbrajQ) turn out to be associated with an end-point contribution 
to the integral at e = 0. 

The basic analysis through section 4 is applicable to a much more general 
case, but the derivation of a relatively simple explicit form of the asymptotic 
Green function, Eq. H96|l . was obtained only for the spherical case. The proce- 
dure can in principle be modified to include cases where the curvature R{t) is 
not piecewise constant. Because the acceleration is not constant, the energy e„ 
in this case is not conserved on the creeping segment. One nevertheless may, 
for instance, express the action Sn{hn) of a classical path with n bounces in 
terms of, say, the maximal "height" /i„ of the classical trajectory. It is not 
difficult to verify that hn is a stationary point of Sn{h) if one allows this height 
to vary while keeping fixed the endpoints and number of "bounces" n. One can 
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then proceed as in section 4 and express the semiclassical approximation to the 
Green function as an integral over h and first sum over the number of bounces 
n before evahiating the integral over h asymptotically. However, unlike the case 
considered in section 4, the summation over n generally cannot be performed 
in closed form if the acceleration is not piecewise constant. One nevertheless 
may sometimes be able to pick out the smallest values of h for which the sum 
over n becomes singular and thus (numerically) obtain the leading asymptotic 
approximation in the general case. 

With a view to applying the method to electromagnetic problems with ideal 
conductors, we considered only Dirichlct boundary conditions, with a change in 
the phase by (—1) at each reflection. Semiclassically, other idealized boundary 
conditions change the phase by a different value. Furthermore, if absorption 
cannot be neglected, the modulus of the reflection coefficient is less than unity. 
Although the physical effect of such changes can be dramatic, since the poles 
are moved, it is a great virtue of the semiclassical approximation that neither 
changes in the boundary conditions nor the inclusion of absorption alter the 
procedure in any fundamental way. 

We briefly comment now on two related papers. Our approach is similar 
in spirit to the one of V.N. BuslaevpS). Buslaev writes the Green function 
G as a "continuum integral" product of Green functions for infinitesimal time 
intervals and uses the saddle point method to obtain the asymptotic form of 
each factor. We do not explicitly use the product form and treat the problem by 
expansion about the classical paths of extremal, but finite, length. Our approach 
appears to be simpler and reduces the diffraction problem to the propagation of 
a particle in a two-dimensional space under the influence of a gravitational-like 
force. However, Buslaev considers a mathematically more general problem than 
we do, that is asymptotic solutions to general parabolic equations. He also does 
not restrict the dimensionality of space to three. Furthermore, our treatment of 
the boundary conditions is very different from that of Budaev. McLaughlin and 
Kellerj23I note that Buslaev's result agrees with Keller's geometrical theory of 
diffraction for the field of a surface diffracted ray. An interesting alternative to 
these methods may be the solution of the wave equation by analytic continuation 
of a random walk proposed in ref . j^S] . 

Acknowledgements: We thank J.B. Keller for bringing references ;23j and|25| 
to our attention and for a number of useful comments. This work was supported 
by the National Science Foundation with Grant PHY-0070525. M.S. greatly 
enjoyed the hospitality of New York University. 

A The Proportionality Factor 

To derive the relation of Eq. 1)28(1 between the semiclassical energy Green func- 
tion G'(X, Y; E) and the Green functions G-y^ defined in Eq. l|77|l that describes 
transverse deviations from a classical path 7c, consider first the general case of 
a medium with a smooth, everywhere differentiable index of refraction 7t.(X, E). 
In this case classical trajectories 7c from X to Y are stationary points of the 
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action and the usual semiclassical formalism applies. In particular, the semi- 
classical energy Greenfunction is given^^ by, 



G(X, Y;E)=Y, \/:D;:exp ( ^^(7^, E) 



27rr 



n 



(112) 



For a given classical path 7c, D^^ is the determinant of a 4 x 4 matrix, 



= det 







axav 


axdE 






dYdE 


dE'-' 



\ 



(113) 



with elements that are 3x3, 3x1, 1x3 and 1x1 matrices. The form of D^^ 
can be simplified by choosing a local coordinate system in the neighborhood 
of the trajectory 7c from X to Y: the coordinate axis for x\\ runs along the 
particular trajectory and the remaining coordinates x = {xrTX±) are transverse 
to the trajectory. Since \/vg at the endpoint of a trajectory is just the change 
in the total time Eq. © associated with an infinitesimal displacement of the 
endpoint in a direction tangent to the trajectory, two of the second variations 
of the action in Eq. I|113() are given by the group velocities at the endpoints, 
namely. 



Vg\Y,E) 



d^Sij,,E) 
dxj dxjdE 

9y\\ dyidE 



In this coordinate system, D^^ becomes [T7| 



Vg{X,E)Vg{Y,E) 



■det 



-d'Si^,,E) 
dxdy 



(114) 



(115) 



proportional to the determinant of a 2 x 2 matrix. Crucial in deriving Eq. H28|l 
is the observation that the matrix of second derivatives of S with respect to 
transverse deviations from the classical path 7c at the endpoints can also be 
obtained by expanding S to quadratic order in the transverse deviations of the 
path 7 from the classical path 7c. We thus have that 



-d'Sij,,E) 
dxdy 



0;r(7c,i?)) 



dadb 



(116) 



On the other hand, if w(X, E) is a smooth function of the coordinate, the clas- 
sical trajectory 7c is stationary and S**"^ is quadratic in the fluctuations. In this 
case, the semiclassical approximation to fj-y^(a, 6;r) is exact. The solution to 
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the equation of motion for a quadratic action S""^ is unique and there is only 
one classical trajectory 7c from a to & in time r that contributes to the two- 
dimensional Greenfunction Q^^{a, 6; r). Van Vleck's formula in two-dimensional 
space then gives 



2Trih 



det(^^5§^)cxp(l5-(a,^;.)^ (117) 



\ \ dadb 



n 



Note that the only classical trajectory from a = to & = in time t{'Jc, E) is the 
trivial one with vanishing transverse deviation everywhere and that therefore 

5(7c,i?)-^;^(5 = 6 = 0;T(7„i?)) . (118) 

Comparing Eq. (|117|) with Eq. (|112|l for r = t{'Jc,E) and a = & = and using 
Eqs. H118|l . (|116|l and l|115|l . one arrives at the relation Eq. JJHl- 

Strictly speaking the validity of Eq. (|28|l has only been shown for a smooth, 
everywhere differentiable phase velocity w(X, S). In this case the classical tra- 
jectories 7c are stationary points of the action and the scmiclassical expressions 
we employed are valid. We argue that Eq. H28() also holds when t;(X, E) is a step 
fimction that vanishes within V, because Eq. H28|) does not explicitly depend on 
the spatial dependence of w(X, E) and involves only the group velocities at the 
endpoints of the classical trajectories. Eq. (|28|l does not depend on the nature 
of the obstacle and the discontinuous spatial dependence of the phase velocity 
of interest can be considered as a limiting case. One similarily can argue that 
because Eq. (|28|) holds for endpoints X and Y connected by classical trajec- 
tories that do not touch the obstacle, the relation should by continuity remain 
valid when one of the endpoints lies in the classical shadow region of the other. 



B Proof that V^iO ^ D^iO for n ^ oo 

It is instructive to see that the diffraction amplitudes defined in Eqs. I|57|l 
and (|86() coincide in the limit of large n for any fixed real value of ^. For 
large n ^ oo the zero's of the Airy function e„ approach e„ defined by Eq. (|83|l : 
in particular, e„ becomes arbitrarily large for n — > oo. 
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Fig. 5; The contour F in the representation for the Airy- function Ai(2) of Eq. 111911 . 
For large values of \t\, F approaches the asymptotes with arg t = ±27r/3 shown as 
dashed lines. 

The Airy function can be represented by the contour integral [TB| 

Ai(z) = — f dte*^"*'/3 , (119) 

where the contour F is sketched in Fig. 5. The expression for I'n(C) of H57|l 
then becomes 



-1 poo p 2 3 /— 



(120) 



For large values of e„ the saddle-point approximation to the integral of Eq. (|120|l 
is accurate. The two saddle points of the integrand arc located at 



p^ = 2^ipt{0, tt = -ipt{0), (121) 



with 



-iO^±yA^+e-^ ■ (122) 



For large values of n, p~ is negative and does not contribute to the saddle 
point approximation of the integral. The path of integration can, however, 
be deformed to pass over the other stationary point at (p^,t^). The integral 
of Eq. (|12()|l in saddle point approximation thus becomes. 



^"(^)~ |Ai'(-e„)|(6„+C^)i/4"^P* 



(123) 



One finally arrives at the expression of Eq. (|86|) for £)„ (^) by using the asymp- 
totic formula for the Airy function. Eq. H47|) implies that for large n, 

en e„, P+(C) PniO and |Ai'(-e„)l - e'/ VV^ , (124) 
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where e„ and Pn(C) are defined in Eqs. (|S5|l and (|S7|l . respectively. Although 
the amplitudes T>n{C) and Dn{C) coincide for large values of n, this does not 
mean that the contribution from the poles of Eq. H85fl , by itself, gives the correct 
asymptotic expansion of the Green function. The asymptotic expansion of the 
Green function is in fact dominated by low-n terms, and we indeed obtained 
asymptotic corrections to the pole contributions in the semiclassical approxima- 
tion. 



C Asymptotics of P{z) 

The function P{z) defined by Eq. (|92|) may be written as a McLaurin series, 

n=l ^ ' 



with coefficients Pn given by 



C, 1 + ie-*«' 



-ad k=l 

(126) 

For given fc, the integral for the contour Csd of Fig. 4 can be performed by 
making the substitution 

? = \/p/ke-"^^ , (127) 
where p is real, and using the definition (|94() . We then have 

Pr^ = /^r(2n/3)f:^ 

= z(2z/3)^/'e-'^^r(2n/3)Li2„/3(-i) . (128) 

Eqs. H125(l and p28|l give the convergent series of Eq. (jnSJ- Note that the 
series in Eq. H93(l uniquely defines P{z) in the whole complex plane but for 
l^l ^ 1, the accurate evaluation of P{z) requires a sizable number of terms and 
is numerically not very efficient. 

For \z\ 3> 1, terms of order in the exponent may be ignored to leading 
order of the steepest descent method and the integral of Eq. (|92(l approaches 
the asymptotic form in Eq. H95|l . 

p{z » 1) ^d^YTT ^ ¥7^z' ■ 

For real x, P{x) is a rather slowly changing function of x that smoothly inter- 
polates between P(0) and P{x 3> 1). For real positive x, P{x) is numerically 
approximated by the rational function 

^(-•-o^TTT^ (130) 

0.7(1 — I) + X 

to within about 10%. 
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